Imaging Neuroscience

Harmonisation is widely used to mitigate site- and scanner-related batch variability in multisite neuroimaging studies and is particularly critical in longitudinal clinical trials, where detection of subtle biological or treatment-related changes depends on reliable measurement across scanners and timepoints. However, the effectiveness of harmonisation in small, heterogeneous clinical datasets remains insufficiently understood, particularly in relation to subject-level variability and consistency across acquisition settings, and its impact on both removal of technical variability and preservation of biological variation in pooled multisite analyses. We systematically evaluated a range of image-based and statistical harmonisation methods using a clinically realistic multisite, multiscanner structural T1-weighted (T1w) MRI test-retest dataset comprising three controlled acquisition scenarios: repeatability, intra-scanner reproducibility and inter-scanner reproducibility. Methods were applied under different batch specifications (site, scanner, or both) and performance was assessed within each scenario and in pooled data using a multi-metric framework capturing both technical and biological variability in volumetric imaging-derived phenotypes (IDPs) relevant to aging and dementia research. Across IDPs, before harmonisation variability was lowest in the repeatability scenario (median variability=0.6 to 2.7%, rank consistency {rho} [≥]0.9), with modest increases under intra-scanner reproducibility (0.5 to 3.2%, {rho}=0.5 to 1.0) and substantially greater variability under inter-scanner reproducibility conditions (1.7 to 19.2%, {rho} =-0.1 to 0.9). These results offer important information to consider for multisite study design, including sample size calculation in clinical trials. Harmonisation performance was strongly context dependent, with clearer benefits emerged in inter-scanner scenarios where both variability reduction and improvements in subject-level consistency were observed. In pooled data, approaches that explicitly modelled site as batch and accounted for repeated-measure structure showed greater consistency across IDPs in batch effect mitigation and more accurately reflected underlying biological variation. Our evaluation metrics enabled disentangling the removal of global batch effect while highlighting residual variability at the phenotype-specific or multivariate levels. These findings demonstrate that harmonisation cannot be treated as a one-size-fits-all solution and must be interpreted relative to the acquisition context, dataset structure, and downstream analytic goals. Multi-metric evaluation under realistic clinical constraints is essential to support reliable and translatable neuroimaging inference by ensuring appropriate correction of batch effects while preserving longitudinal biological signals and sensitivity to clinically meaningful change in multisite studies.

The Adolescent Brain Cognitive Development (ABCD) Study is the largest U.S.-based neuroimaging initiative of adolescent brain maturation. Diffusion MRI (dMRI) provides unique insights into white matter organization, yet applying advanced processing pipelines and managing technical variability across scanning environments remains challenging at scale. To address these issues, we present ABCD-BIDS Community Collection (ABCC) release 3.1.0, including a curated resource of more than 24,000 fully processed ABCD dMRI datasets. ABCC provides fully processed images, nuanced image quality metrics, advanced microstructural measures, and person-specific bundle tractography. Evaluating these rich data revealed that measures of diffusion restriction and non-Gaussianity--in particular the intracellular volume fraction from NODDI and return-to-origin probability from MAP-MRI--were highly sensitive to neurodevelopment and robust to variation in image quality. Additionally, harmonization of microstructural features markedly improved the cross-vendor generalizability of developmental effects. Together, ABCC accelerates reproducible, rigorous research on adolescent white matter development.

Multi-site diffusion MRI studies face scanner-induced variability that can obscure biological signal. Harmonization methods such as ComBat have been developed to address this, but have been evaluated primarily on diffusion scalar metrics. Whether scanner reproducibility differs across fundamentally distinct tract-derived representations has not been systematically compared. Here, we compared the reproducibility of three metric families (diffusion, shape, and connectivity) across 36 association tracts using the MASiVar dataset (5 subjects, 4 scanners, 27 sessions). We assessed intraclass correlation coefficients (ICC) and multivariate subject discrimination at baseline, under dimensionality reduction, and after ComBat harmonization. At baseline, shape metrics showed the highest reproducibility (median ICC 0.69), followed by connectivity (0.49) and diffusion (0.34). Shape and connectivity achieved comparable subject discrimination (both 1.75), significantly exceeding diffusion (1.23). ComBat harmonization improved all families but harmonized diffusion (0.58) remained below unharmonized shape (0.69), indicating that metric family selection remains consequential even after harmonization. Under low-dimensional representation, connectivity showed the largest gains (ICC 0.86, subject discrimination 3.0), exceeding other families at any dimensionality. Analysis of principal component loadings identified a small number of cortical regions per tract (median 6) that capture 95% of the reproducible connectivity signal, providing a per-tract reference for selecting the most informative regions in future multi-site studies. These findings indicate that the choice of which tract-derived metrics to analyze in multi-site studies deserves at least as much consideration as how to harmonize them.

Subcallosal cingulate cortex (SCC) deep brain stimulation (DBS) can provide relief for individuals with Treatment Resistant Depression (TRD), but ongoing clinical management remains challenging due to nonspecific symptom fluctuations that can obscure core depression recovery on standard rating scales. Objective, stable biomarkers that selectively track the therapeutic effects of SCC DBS are therefore essential for developing principled decision support systems to guide stimulation adjustments. Recent bidirectional DBS systems enable chronic recording of local field potentials (LFPs) and prior work using the Activa PC+S device identified an electrophysiological signature of stable clinical recovery. However, translation to practical clinical deployment requires demonstrating that this biomarker is robustly generalizable, specific to the impact of the DBS therapy, and deployable in real-world recording contexts. To address this need, we developed an at-home SCC LFP data collection platform (built on the Medtronic Summit RC+S system) enabling at home data collection for a new cohort of ten SCC DBS participants with TRD (ClinicalTrials.gov identifier NCT04106466). Using longitudinal LFP recordings collected from this system, we report findings demonstrating that the previously reported biomarker of stable recovery generalizes across subject cohorts and devices, is robust to common potential confounds (including time of day and stimulation status), and shows symptom specificity, sensitivity and stability necessary to support clinical decision making. Across both cohorts, biomarker changes show relationships to pre-DBS white matter structure and network function measured using diffusion MRI and resting-state functional MRI (rsFMRI). These findings replicating and extending previous findings support the biomarkers utility as a foundation for scalable, electrophysiology-informed decision support in SCC DBS.

Electroencephalography (EEG) preprocessing is a critical yet time-consuming step that often relies on expert-driven, semi-automatic pipelines, limiting scalability and reproducibility across large datasets. In this work, we present sEEGnal, a fully automated and modular pipeline for EEG preprocessing designed to produce outputs comparable to expert-driven analyses while ensuring consistency and computational efficiency. The pipeline integrates three main modules: data standardization following the EEG extension of the Brain Imaging Data Structure (BIDS), bad channel detection, and artifact identification, combining physiologically grounded criteria with independent component analysis and ICLabel-based classification. Performance was evaluated against manual preprocessing performed by EEG experts at two complementary levels: preprocessing metadata (bad channels, artifact duration, and rejected components) and EEG-derived measures. In addition, test-retest analyses were conducted to assess the stability of the pipeline across repeated recordings. Results show that sEEGnal achieves performance comparable to expert-driven preprocessing while preserving key neurophysiological features. Furthermore, the pipeline demonstrates reduced variability and increased consistency compared to human experts. These findings support sEEGnal as a robust and scalable solution for automated EEG preprocessing in both research and large-scale applications. HighlightsFully automated and modular EEG preprocessing pipeline. Benchmarked against expert-driven preprocessing. Comparable performance in metadata and EEG-derived measures. Demonstrates stable performance in test-retest recordings. BIDS-based framework for reproducible EEG data handling.

ObjectiveSocial difficulties are transdiagnostic in childhood, but their heterogeneity is poorly characterised and rarely treated as a primary neurodevelopmental phenotype. This matters because childhood and adolescence are sensitive periods for peer relationships and brain development. We used data-driven modelling and non-linear mapping to derive social profiles and test their clinical, cognitive, and neural correlates. MethodsParticipants were 992 children aged 5-18 years from CALM (Mage = 9.6). Social items from the SDQ, CCC-2, and Conners-3 were modelled using a regularised partial correlation network to derive core social dimensions. A self-organising map captured graded social profiles. Simulated archetypes, SVM-based island identification, and permutation testing defined profile regions and centroid-distance scores. Profiles were related to referral, diagnosis, cognition, BRIEF indices, and T1-derived MIND network structure in an MRI subsample (n = 431). ResultsWe identified four profiles: social engagement, friendship difficulties, social withdrawal, and peer victimisation. Profile expression tracked variation in referral and diagnostic pathways. Social withdrawal showed the clearest disadvantage across cognitive domains, whereas social engagement was associated with fewer executive function difficulties across BRIEF indices. MIND strength components covaried with profile expression (a significant PLS latent variable, p = 0.02), with covariance strongest for social withdrawal and peer victimisation. ConclusionsChildhood social functioning organises graded signatures that relate to clinically relevant pathways, cognitive and executive outcomes, and brain structure. Profiling social signatures provides a scalable framework for identifying social need beyond diagnostic categories, motivating studies to test directionality and improve developmental outcomes.

ObjectivePathological beta-band oscillations (13 to 30 Hz) in the subthalamic nucleus (STN) are a hallmark of Parkinsons disease and a primary target for deep brain stimulation therapy, yet the specific pattern of synaptic reorganization that drives their emergence remains incompletely understood. We developed a GPU-accelerated computational framework to systematically investigate combinations of synaptic changes across basal ganglia pathways that produce Parkinsonian beta oscillations while satisfying literature-based electrophysiology constraints. ApproachWe implemented a biophysically detailed spiking network model of the STN, external globus pallidus (GPe), and internal globus pallidus (GPi) in JAX (a high-performance numerical computing Python library), achieving a 490-fold speedup over conventional CPU-based simulation. Using the Covariance Matrix Adaptation Evolution Strategy (CMA-ES) we optimized 10 network parameters across two stages: first establishing a healthy baseline matching primate electrophysiology data, then searching within biologically motivated bounds for synaptic modifications that reproduce Parkinsonian firing rates and beta power. Fixed in-degree connectivity ensured optimized parameters produced scale-invariant dynamics from 450 to 45000 neurons. All simulations ran on a single cloud GPU instance at 84 cents per hour. Main ResultsThe optimizer converged on a coordinated pattern of synaptic reorganization dominated by asymmetric changes within the STN-GPe reciprocal loop: STN to GPe excitation increased 2.21-fold while GPe to STN inhibition collapsed to 0.11-fold of its healthy value. STN to GPi and GPe to GPi pathways changed minimally (1.06-fold and 1.45-fold respectively). This configuration transformed asynchronous firing (beta: 0.4 percent of spectral power) into synchronized bursting with prominent beta oscillations (49.4 percent), with firing rate changes matching experimental observations. Network dynamics were invariant across a 100-fold range of network sizes (firing rate deviation less than 2.4 Hz; all metrics p less than 0.001 across 10 random seeds at 45000 neurons). We implemented a simplified deep brain stimulation model for validation purposes, which achieved complete beta suppression (49.4 percent to 0.0 percent) and restored GPi output to healthy levels. SignificanceThese results suggest that pathological beta oscillations emerge from a specific pattern of synaptic reorganization, namely the reduction of GPe inhibitory feedback to STN. The GPU-accelerated optimization framework, running on commodity cloud infrastructure, demonstrates an accessible platform for parameter exploration in neural circuit models and a foundation for generating synthetic training data for adaptive deep brain stimulation algorithms.

Understanding neural correlates of brain function in neuroscience now largely involves detecting and analyzing transient signals from fluorescent sensors. Imaging technologies such as confocal and two-photon microscopy, along with onboard miniscopes, enable the visualization of neural activities and capture dynamic signals both ex vivo and in vivo. This includes monitoring Ca2+ transients via the expression of genetically encoded sensors such as GCaMP in specific brain cells. Additionally, the advent of GPCR-based neurotransmitter sensors allows for imaging the release of neurotransmitters including glutamate and GABA, as well as neuromodulators such as dopamine or noradrenaline. These approaches however generate large, high-dimensional, spatiotemporally complex datasets, presenting significant challenges for signal detection and analysis. To overcome these challenges, we developed a versatile pipeline of Dynamic Extraction and Tracking of Emitted Cellular Transients (DETECT), which combines background denoising, object segmentation, and multi-object tracking. Our user-friendly, Python-based GUI offers a low-resource platform for efficient data analysis. Validated across various imaging modalities and biological models, DETECT provides a robust and comprehensive solution for analyzing complex imaging datasets in neuroscience research.

Brain network analysis has emerged as a critical framework for understanding the complex organization and function of the human brain, underpinning insights into cognition, behavior, and neuropsychiatric conditions. Central to this approach is the parcellation of the brain into discrete regions, which simplifies high-dimensional connectome data and facilitates the investigation of network architectures. However, the proliferation of brain parcellation schemes introduces significant challenges: different parcellations often yield varying network sizes and measures, complicating cross-study comparisons and the reproducibility of findings. Moreover, most connectome construction pipelines are rigid, typically outputting connectivity matrices from only one or a few parcellation schemes, which limits flexibility. In this paper, we address these issues by introducing BridgeBP, a novel toolbox designed to bridge brain parcellations by leveraging continuous brain connectivity concepts. BridgeBP transforms structural connectivity matrices derived from one parcellation scheme into matrices corresponding to more than 40 alternative schemes, standardizing analyses and enhancing the robustness of network studies. Through extensive evaluations, we demonstrate that BridgeBP enables consistent network comparisons across diverse parcellation frameworks, paving the way for more reproducible and generalizable insights in brain connectome research.

Cognitive reappraisal, the deliberate reinterpretation of emotional events, is widely considered an effective emotion regulation strategy, and modulation of the late positive potential (LPP) during negative affect reduction has become the primary electrophysiological evidence for volitional emotional control. Experimental instructions, however, impose dual-task demands that free viewing does not, confounding reappraisal with cognitive load. By including instructions to increase emotional responses to pictures ("enhance") as well as instructions to decrease ("suppress"), different predictions are generated. If the LPP reflects regulation, then, compared to free viewing, suppress instructions should decrease LPP amplitude, and enhance instructions should increase LPP amplitude. If modulation instead reflects cognitive load, both instructions should reduce the LPP, as both impose an additional cognitive task. In a sample of 107 participants, evaluative ratings confirmed that regulation instructions modulated reported emotional intensity in the expected directions (Enhance > View > Suppress), but that both enhance and suppress instructions reduced LPP amplitude compared to free viewing, with Bayesian model comparisons providing strong evidence against direction-specific regulation and in favor of cognitive load. Whole-scalp multivariate pattern analysis confirmed that no instruction-related neural signal exists at any scalp location or latency within the first second after stimulus onset. These data indicate that LPP modulation following both instruction types reflects dual-task cognitive load rather than volitional emotional control. Significance StatementCognitive reappraisal is considered the gold standard of emotion regulation, and reduced late positive potential (LPP) amplitude during negative emotion suppression is the primary neural evidence that humans can voluntarily control emotional responses. The current data are inconsistent with this regulatory account and instead support a cognitive load interpretation. Whether instructed to enhance or suppress emotional responses, LPP amplitude was reduced in both conditions relative to free viewing, consistent with attentional resource competition rather than directional regulatory control. The same participants reported successfully regulating emotional experience in opposite directions, producing a clear dissociation between neural and behavioral measures. These findings challenge a basic tenet of emotional regulation and raise questions concerning LPP modulation as a biomarker of regulatory capacity.

Hidden (semi-) Markov Models (HsMMs) are increasingly being used to segment neurophysiological signals into sequences of latent cognitive processes. The idea: different processes will leave distinct traces in trial-level recordings of (multivariate) neuro-physiological signals. Markov models, equipped with an emission model of these traces and a latent process model describing the progression through the different latent processes involved in a task, can then be used to infer the most likely process for any time-point and trial. However, the currently used HsMMs remain limited in two important ways. First, they cannot account for subject-level heterogeneity in the latent and emission process. Instead, a single group-level model is assumed to explain the entire data. Second, they cannot account for the potentially non-linear effects of experimental covariates on the latent and emission process. To address these problems, we present a modeling framework in which the HsMM parameters of the emission and latent process are replaced with mixed additive models, including smooth functions of experimental covariates and random effects. We derive all necessary quantities for empirical Bayes and fully Bayesian inference for all parameters and provide a Python implementation of all estimation algorithms. To demonstrate the advantages offered by this framework, we apply such a multi-level model to an existing lexical decision dataset. We show that, even in such a simple task, not all subjects rely on the same processes equally and that at least two semi-Markov states, previously believed to reflect distinct processes, might actually relate to the same cognitive process.

Accurate segmentation of nasal and paranasal sinus structures from CT scans is critical for surgical planning and treatment evaluation in rhinology. However, the complex anatomical topology and thin-wall boundaries of these structures pose significant challenges for automated segmentation methods. We propose AFS-DSN (Adaptive Frequency-Spatial Dual-Stream Network), a novel deep learning architecture that integrates multi-scale wavelet decomposition with spatial feature learning for binary segmentation of the nasal cavity complex. Our method employs a dual-stream encoder with frequency branch utilizing three wavelet scales (db1, db2, db4) to capture 24 frequency sub-bands, enabling enhanced boundary detection in anatomically challenging regions. Cross-domain attention and adaptive routing mechanisms dynamically fuse spatial and frequency features based on local tissue characteristics. We formulate the task as binary segmentation where all five anatomical structures (maxillary sinus, sphenoid sinus, ethmoid sinus, frontal sinus, and nasal cavity) are treated as a unified foreground region against the background, prioritizing clinical boundary detection over individual structure differentiation. Evaluated on the NasalSeg dataset (130 CT volumes) with a 70/15/15 train/validation/test split, AFS-DSN achieves 94.34% {+/-} 2.30% overall Dice coefficient with statistically significant improvements in thin-wall regions (91.34% vs. 90.57% baseline, p=0.004) and statistically significant improvement in Surface Dice at 1mm tolerance (0.874 vs. 0.868 baseline, p=0.010), demonstrating enhanced boundary precision while maintaining sub-second inference time, making the method suitable for surgical planning applications where sub-millimeter accuracy is clinically relevant. To address concerns regarding model complexity, we further introduce AFS-DSN-Lite, a parameter-efficient variant (27.41M parameters) that achieves comparable performance (94.37% Dice) through depthwise separable convolutions, and validate robustness via 3-fold cross-validation (mean Dice: 94.59% {+/-} 0.31%).

Alzheimer's disease (AD) is increasingly recognized to have systemic physiological correlates alongside central neurodegeneration. Here, we explored brain-organ network (BON) connectivity in AD (n=28) and healthy controls (n=23) using time-resolved quasi-dynamic analysis of plateau-phase total-body 18F-tau-PET. We found that AD-related pathophysiology was linked not only to cerebral tau aggregation, but also to altered signal synchronization across the brain-organ network, despite comparable body tracer distribution. Network topology analyses revealed the occipitotemporal cortex and the spinal cord as key nodes in this altered systemic network. Furthermore, exploratory mediation analyses demonstrated that BON dysregulation is cross-sectionally linked to cognitive deficits, with statistical associations observed for both cortical tau burden and imaging markers of impaired glymphatic clearance. This total-body PET study provides first-ever direct evidence repositioning AD as a multi-organ disorganization disease. These findings provide a novel framework for investigating brain-body interactions and systemic vulnerabilities in neurodegenerative disorders.

Chronotype reflects individual circadian preference for timing of sleep, wakefulness, and peak performance and has been linked to variability in prefrontal cognitive function across the day. Whether chronotype independently relates to dual-task gait cost (DTC) and whether this relationship differs by cognitive task domain is unclear. Sixty-nine healthy young adults (37 female; mean age 21.3 years) completed the Morningness-Eveningness Questionnaire (MEQ). Spatiotemporal gait parameters were recorded with three-dimensional motion capture during single-task walking and three dual-task conditions: backward word spelling (5LWB; phonological), serial subtraction by seven (SS7; arithmetic), and reverse month recitation (RMR; sequential). DTC was calculated for eight gait parameters. Condition differences were assessed with nonparametric tests and post-hoc comparisons. Multiple linear regression, adjusting for age, sex, BMI, and baseline gait velocity, tested the independent association between MEQ score and mean velocity DTC; exploratory Spearman correlations examined other parameters. SS7 produced the largest mean velocity DTC (-12.76%), significantly greater than 5LWB (-7.95%; p = 0.002) and RMR (-9.57%; p = 0.021). MEQ score independently predicted mean velocity DTC in 5LWB ({beta} = -0.51, p < 0.001, R{superscript 2} = 0.269) and RMR ({beta} = -0.55, p = 0.004, R{superscript 2} = 0.222), indicating greater morningness associated with better gait-speed preservation under cognitive load; the SS7 association was not significant ({beta} = -0.33, p = 0.071). Exploratory correlations showed MEQ-DTC associations across 7/8 parameters in 5LWB, 4/8 in RMR, and 3/8 in SS7. Chronotype is independently associated with dual-task gait cost in a task-domain-specific manner, with stronger effects for phonological and sequential tasks than for arithmetic processing. The SS7 condition yielded the largest interference but weakest chronotype modulation, suggesting arithmetic dual-task disruption may be less sensitive to circadian arousal. Fixed testing time and cross-sectional design warrant within-subject, multi-timepoint studies to confirm chronotype effects separate from time-of-day confounds.

This research demonstrates that the trans-aqueduct approach is a feasible, minimally invasive access pathway to the third ventricle, offering a potential route to the deep brain for therapeutic technologies. Further pre-clinical investigation is required to thoroughly evaluate physiological tolerance, trauma risk, and the long-term implications of intraventricular implantation. The third ventricle is a high-value site for neuromodulation due to its proximity to deep-brain targets, including the subthalamic nucleus (STN) and globus pallidus internus (GPi). This study defined the anatomical pathway; and evaluated the technical feasibility of retrograde access to the third ventricle via the cerebral aqueduct using minimally invasive interventional techniques. Evaluation was conducted in three phases using human MRI datasets (n=16; mean age 48.4 years) and cadaveric specimens (n=6; mean age 88.2 years). Phase 1 involved morphometric MRI analysis of the aqueduct and ventricles. Phase 2 tested trans-aqueduct access on cadaver specimens via fluoroscopically guided guidewires and catheters. Phase 3 utilized direct anatomical dissections on cadaver specimens (n=3) to morphometrically measure the third ventricular cavity and its relationship to deep-brain nuclei. Measurements across the sample groups showed a mean aqueduct diameter of 1.6 mm (SD=0.14). Third ventricle dimensions averaged 27.6 mm (ventral-dorsal), 19.9 mm (caudal-cranial), and 5.7 mm (lateral). Successful access to the third ventricle was achieved in 83% (5/6) of cadaveric specimens. The optimal technical configuration utilized a 0.018'' angled-tip guidewire and 5-6 Fr catheters; the aqueduct accommodated diameters up to 2.0 mm with minimal resistance. The STN and GPi were localized within 5-20 mm of the ventricular volumetric centroid. The trans-aqueduct approach is a technically feasible, minimally invasive pathway for accessing the third ventricle. This route offers a potential alternative for the delivery of therapeutic neurotechnologies. Further research is required to assess physiological tolerance, trauma risk, and the long-term safety of intraventricular implantation.

Extracellular vesicles (EVs) mediate intercellular transfer of lipids, proteins, and nucleic acids between nearly all cell types. We previously showed that astrocyte-derived EVs modulate neuronal mitochondria in vitro, but whether endogenous astrocytic EVs are trafficked to neuronal mitochondria in vivo remained unknown. To address this, we generated an EV reporter mouse, Aldh1l1-Cre; CD9-tGFPfl/fl, in which astrocyte-secreted EVs are labeled with a CD9-turboGFP fusion protein (CD9-tGFP). Astrocyte-specific expression of CD9-tGFP was verified in brain tissue and isolated EVs, comprising 13.2 {+/-} 1.6% of total brain EVs. In primary glial cultures, CD9-tGFP was restricted to astrocytes, localizing to vesicular compartments and cell protrusions (filopodia and cilia), with 89.3 {+/-} 2.2% of astrocyte-derived EVs carrying the label. These EVs were enriched with the sphingolipid ceramide, consistent with its co-distribution with CD9-tGFP in astrocytic cell protrusions. In the cortex, hippocampus, and cerebellum, CD9-tGFP was predominantly detected in astrocytic processes co-labeled with GLAST1 and GFAP, forming contacts with laminin-positive capillaries and parvalbumin-positive neurons. CD9-tGFP-labeled EVs were detected inside capillaries and neurons, and super-resolution STED microscopy revealed partial overlap with neuronal mitochondria. Live-cell spinning disk confocal imaging and AI-assisted proximity analysis confirmed uptake of CD9-tGFP EVs by neuronal cells and trafficking of their cargo to mitochondria in vitro. Biochemical isolation of synaptic and non-synaptic mitochondria confirmed EV-derived cargo on mitochondria in vivo, with 3-fold higher association of CD9-tGFP with synaptic than non-synaptic mitochondria. Together, these findings validate the Aldh1l1-Cre; CD9-tGFPfl/fl reporter mouse as a powerful tool for tracking astrocyte-derived EVs in vivo and provide direct evidence that their cargo is preferentially trafficked to synaptic mitochondria. Graphical AbstractAstrocyte-derived extracellular vesicles target neuronal mitochondria in vivo O_FIG O_LINKSMALLFIG WIDTH=156 HEIGHT=200 SRC="FIGDIR/small/718987v1_ufig1.gif" ALT="Figure 1"> View larger version (33K): org.highwire.dtl.DTLVardef@174d92aorg.highwire.dtl.DTLVardef@5d8248org.highwire.dtl.DTLVardef@114483borg.highwire.dtl.DTLVardef@924d55_HPS_FORMAT_FIGEXP M_FIG C_FIG

Early behavioral and temperamental differences are important indicators of later socioemotional development and psychopathology risk, yet their neural bases near birth remain incompletely understood. Using resting-state fMRI data from the Developing Human Connectome Project, we examined whether neonatal functional connectivity predicts 18-month behavioral and temperament outcomes in 397 infants (277 term-born, 120 preterm-born). Outcomes were assessed using the Child Behavior Checklist (CBCL) and the Early Childhood Behavior Questionnaire (ECBQ). We applied a stability-driven, ROI-constrained connectome-based predictive modeling framework to identify robust whole-brain connectivity features associated with later externalizing, internalizing, surgency, negative affect, and effortful control. Significant predictive models were observed for multiple outcomes across the whole cohort as well as within term-born and preterm-born groups, with clear differences in predictive architecture between cohorts. Across analyses, prefrontal and temporoparietal systems were repeatedly implicated, alongside medial temporal, fusiform, parahippocampal, and orbitofrontal-related regions. These findings indicate that large-scale neonatal functional organization is meaningfully related to later socioemotional and behavioral variation, and that preterm birth is associated with partly distinct predictive connectivity patterns.

Salient sensory stimuli are known to evoke neural activations across distributed brain regions. However, the temporal dynamics of these responses over sub-second timescales remain poorly understood, in part due to limitations in the temporal resolution of non-invasive neuroimaging methods. We examined the spatiotemporal dynamics of neural activations evoked by salient sensory stimuli (rare sounds) using 1,194 widely distributed intracranial electrodes in 5 neurosurgical patients. Salient stimuli preferentially activated 263 of 1,194 electrodes (22%), with responses segregating into two largely distinct spatiotemporal patterns: (1) phasic activation in sensorimotor regions, and (2) sustained activation within the salience network. Cross-correlation analysis revealed that phasic sensorimotor activation preceded sustained salience network activation on a trial-by-trial basis. These findings support an updated view of salience processing in the human brain, revealing that salient stimuli evoke two sequential stages of neural activation--phasic sensorimotor responses followed by sustained salience network activity--rather than simultaneous widespread activation.

Background: Transcranial temporal interference stimulation (tTIS) is an emerging noninvasive neuromodulation approach that enables focal, frequency-specific modulation of deep brain regions, offering a novel method for investigating therapeutic mechanisms underlying brain and mental health disorders. Pain is a key target because it is a feature of multiple disorders and is increasingly understood to depend on brain circuits. Here, we tested the effects of tTIS on bilateral evoked pain, capitalizing on converging evidence from human and animal studies indicating that the primary motor cortex (M1) contains body-wide inter-effector regions and has descending projections to regions implicated in nociceptive, motivational, and autonomic processing, making it a key cortical target for pain modulation. Methods: We conducted a pre-registered, triple-blind, randomized crossover study (N = 32, 160 study sessions), investigating frequency-dependent effects of tTIS applied to the left M1 on experimentally evoked thermal pain in healthy adults. We tested four stimulation frequencies (10 Hz, 20 Hz, 70 Hz, and sham) on separate days (>10,000 pain trials total). Noxious heat was applied to both the right and left forearms using individually calibrated temperatures both pre- and post-stimulation. Results: Active tTIS produced significant analgesia at all stimulation frequencies (10 Hz, 20 Hz, and 70 Hz) relative to sham (Cohens d = 0.46-0.82, all p < 0.05). 10 Hz produced the greatest reduction (d = 0.82), and both 10 Hz and 20 Hz produced more analgesia than 70 Hz (d = 0.44 and 0.38, respectively; p < 0.05). Stimulation-related sensations were equivalent across frequencies, and participants were blind to condition. Pain reductions remained stable over a [~]40-min post-stimulation period and were bilateral, consistent with stimulation of body-wide inter-effector regions. Conclusions: These results provide the first evidence that tTIS can reliably reduce experimental pain perception in humans in a frequency-dependent manner, providing a foundation for noninvasive pain modulation with tTIS.

Regional volumetric assessment of perinatal brain development is currently limited by the lack of consistent high quality multi-regional segmentation methods applicable to both fetal and neonatal MRI. We present Multi-BOUNTI, a deep learning pipeline for automated multi-lobe segmentation of fetal and neonatal T2w brain MRI. The method is based on a dedicated 43-label parcellation protocol and a 3D Attention U-Net trained on brain MRI datasets of subjects spanning 21-44 weeks gestational/postmenstrual age. The pipeline integrates preprocessing, segmentation and volumetric analysis, and was evaluated on independent datasets, demonstrating fast (< 10 min/case) and accurate performance with high agreement to manually refined labels. We demonstrate the application of the framework with 267 fetal and 593 neonatal MRI datasets from the developing Human Connectome Project without reported clinically significant brain anomalies to derive normative volumetric growth models across 21-44 weeks GA/PMA. These models were used to characterise developmental trajectories, assess differences between fetal and preterm neonatal cohorts, and analyse longitudinal changes. The resulting normative models were integrated into an automated reporting framework enabling subject-specific volumetric assessment via centiles and z-scores. Multi-BOUNTI provides a unified and scalable approach for perinatal brain segmentation and volumetry, supporting large-scale studies and facilitating future clinical translation. The full pipeline is publicly available at https://github.com/SVRTK/perinatal-brain-mri-analysis.